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INTRODUCTION 


A  series  of  efforts  has  been  made  to  estimate  stress  Intensity  factors  of 
radial  cracks  In  a  cylindrical  pressure  vessel.  In  order  to  facilitate  the 
finite  element  computations  of  such  a  cracked  cylinder,  a  collapsed  12-node 
triangular  element  with  nodes  shifted  to  1/9  and  4/9  locations  was  developed 
in  Reference  l.  The  use  of  these  elements  around  a  crack  tip  has  helped  In 
computing  stress  intensity  factors  for  an  array  of  uniformly  distributed 
radial  cracks  of  equal  depth  emanating  from  the  Inner  surface  of  a 
thick-walled  tube.  Our  finite  element  results,  reported  In  Reference  2,  are 
In  good  agreement  with  results  obtained  by  other  methods  (refs  3-5).  To 
reduce  stress  intensity  factors  at  these  cracks  so  the  fatigue  life  of  the 
pressure  vessel  may  be  prolonged,  it  becomes  a  common  practice  to  use  the 
autofrettage  process.  This  process  produces  a  residual  stress  in  a  cylinder 
following  plastic  deformation  and  elastic  unloading.  The  residual  stress  is 
compressive  near  the  bore  and  is  tensile  near  the  outer  surface.  While  it  Is 
useful  to  retard  the  growth  of  cracks  near  the  bore  where  Che  residual 
compressive  stress  reduces  the  maximum  tensile  stress  level  due  to  a  bore 
pressure,  it  may  cause  problems  at  cracks  near  the  outer  surface  where  the 
final  stress  becomes  the  sum  of  residual  tensile  stress  and  the  tensile  stress 
due  to  a  pressure  loading  applied  at  the  bore.  There  must  be  an  optimal 
degree  of  autofrettage  when  failures  from  both  inner  and  outer  cracks  ire 
equally  likely.  The  computations  of  stress  Intensity  factors  at  Inner  and 
outer  cracks  are  more  involved  if  a  cylinder  undergoes  a  partial  autof rettage. 

References  are  listed  at  the  end  of  this  report. 
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One  of  Che  difficulties  Is  the  discrepancy  among  predictions  on  residual 
stress  distribution  hv  different  Investigators  based  on  different  assumptions. 
References  6  through  13  represent  a  partial  list  of  works  on  this  subject. 
Another  difficulty  Is  the  lack  of  an  existing  general  purpose  finite  element 
computer  code  which  has  the  feature  of  treating  an  arbitrary  initial  stress 
distribution.  To  develop  a  method  of  computing  stress  intensity  factors  in 
this  situation,  the  closed  form  expressions  for  residual  stresses  for  an 
incompressible,  elast Ic-ideally  plastic  material  (ref  7)  are  chosen.  In 
addition,  an  equivalent  thermal  load  was  found  in  Reference  14  which  can  be 
prescribed  at  the  nodes  of  the  existing  computer  codes  NIASTRAN  and  APES  (ref 
15)  to  avoid  the  need  of  modification  of  computer  codes  In  order  to  handle 
Initial  stresses.  A  summary  of  these  methods  and  some  results  in  stress 
intensity  computations  for  cracks  In  a  partially  autof rettaged  cylinder  Is 
given  In  Reference  lb.  Other  efforts  on  this  subject  are  given  In  References 
17  and  13.  Further  extension  of  the  method  developed  for  radial  cracks  In 
cylindrical  pressure  vessels  of  elast lc-oerl:ectly  plastic  'materials  to  similar 
cracks  of  straln-har  leutng  materials  Is  published  In  Reference  IP. 

Thus  fir  the  stress  Intensity  factors  have  been  obtained  for  uniformly 
distributed  radial  cracks  of  equal  lepth.  The  weakest  conf Iguratton  was 
found  to  be  two  diametrically  opposite  cracks  In  non-aut.of rettaged  cylinders. 
For  nutof  retr.  iged  cylinders,  the  weakest  configuration  Is  deoendent  on  the 
degree  of  autofrettage  and  the  ratio  of  bore  pressure  t>  yield  strength.  It 
remains  true  la  general  that  two  diametrically  opposite  cracks  make  a  cylinder 
the  least  resistant  to  fracture.  however,  from  the  actual  observations  of 


fracture  of  thick-walled  tubes,  it  Is  a  single  major  crack  which  has  caused 


the  failure  most  frequently.  The  objective  of  this  study  Is  to  compute  the 
change  In  stress  Intensity  factors  at  different  crack  tips  when  they  grow  from 
an  equal  depth  configuration  to  an  unequal  depth  configuration.  Based  on  some 
finite  element  computations  of  stress  Intensity  factors,  a  numerical  method  is 
needed  to  estimate  stress  Intensity  factors  at  each  crack  tip  of  a  cracked 
tube  when  radial  cracks  are  unequal  depths.  Results  in  this  report  are 
limited  to  non-autof rettaged ,  pressurized  tubes.  However,  it  Is  not  difficult 
to  extend  the  method  to  partially  autof rettaged  tubes.  Results  will  appear  in 
a  separate  report. 

FINITE  ELEMENT  RESULTS  OF  STRESS  INTENSITY 

The  12-node  quadrilateral  isoparametric  elements  were  used  in  this  study. 

The  shape  functions  and  the  mathematical  formulation  of  the  element  stiffness 
matrix  are  given  in  Reference  l.  The  crack  tip  elements  are  formed  by 
collapsing  the  quadrilateral  elements  into  triangular  elements  around  the 
crack  tip.  On  each  side  of  a  triangle  emitting  from  the  crack  tip,  the 
locations  of  Intermediate  nodes  are  shifted  from  their  usual  1/3  and  2/3  of 
the  length  of  the  side  to  1/9  and  4/9  of  the  length  measured  from  the  tip. 

This  simple  technique  gives  the  required  singularity  of  the  strain  field  at 
the  crack  tip  (ref  1).  The  stress  intensity  factors  can  be  computed  quite 
accurately  from  certain  nodal  displacements  (ref  1).  A  general  purpose 
computer  program  having  Isoparametric  elements  such  as  NASTRAN  may  be  used  for 
our  study.  However,  the  computer  code  APES  is  chosen  due  to  many  convenient 
features  of  the  program. 
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Us lnf»  similar  notations  as  In  deference  2,  n  denotes  the  number  of  radial 
cracks,  W  Is  the  wall  ratio,  R2/R1,  where  and  R2  are  Inner  and  outer  radii 
of  the  ring,  respectively.  The  thickness  of  the  wall,  t  =  R2  -  Rj ,  is  used  to 
normalize  the  crack  depth  a.  The  dimensionless  crack  depth  a/t  is  denoted  by 
c.  A.  subscript  l,  l  =  l,2,...n.  Is  used  to  name  a  crack.  Let  the  Cartesian 
coordinates  be  chosen  such  that  the  origin  is  at  the  center  of  the  ring  and  a 
radial  crack  lies  on  the  positive  x-axls.  This  crack  Is  assigned  the  number 
one.  The  crack  number  increases  counterclockwise  from  the  x-axls.  In  this 
study  our  finite  element  computations  are  limited  to  Si  *  2  and  n  =  2,  3  and  4 
as  shown  la  Figure  1.  Starting  from  equal  depth  crack  configuration,  we 
proceed  to  unequal  depth  conf Igurat ions  by  a  systematic  increment  of  crack 
depth  of  one,  two,  or  three  cracks.  For  a  given  n,  let  c  and  Ke  be  the  crar 

depth  and  stress  Intensity  factor,  respectively,  when  all  cracks  are  of  the 

same  depth.  These  quantities  are  used  to  normalize  the  crack  depth  cj  and  the 
mode  1  stress  intensity  factor  '<£  for  the  i-th  crack  when  cracks  become 
unequal  depths.  Vow  notations  e  and  V  are  introduc’d  for  these  dimensionless 
quantities 

v;i  -  c.L/c  ,  Mt  -  Kt/Ke  i  =  l,2,...n  (1) 

It  Is  true  that  some  cracks  are  under  mixed  mode  I  and  mode  II  conditions  when 

unequal-depth-crack  conf igurat tons  are  considered,  however,  in  this  study  the 
effect  of  the  shear  mode  Is  assumed  to  be  negligible  and  only  the  dominant 
mode  I  stress  intensity  factors  are  considered. 

The  idealizations  used  in  this  study  for  a  ring  with  two,  three,  and  four 
cracks  are  shown  I  1  figures  2,  3,  and  4,  respectively.  Only  the  upper  half  of 
a  ring  Is  shown  due  to  symmetry  consideration  In  finite  element  computations. 
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Equilateral  triangles  are  used  around  a  crack  tip  and  these  collapsed  singular 
elements  are  surrounded  by  a  layer  of  quadrilateral  elements,  detail  A  of 
Figure  2  and  detail  8  of  Figure  3.  The  finite  element  results  of  Ke  using  the 
refined  meshes  are  shown  in  Figure  5  and  are  very  close  to  results  previously 
obtained  in  Reference  2  using  slightly  different  meshes  around  a  crack  tip. 

For  two  diametrically  opposed  radial  cracks,  stress  intensity  factors  are 
computed  for  various  lengths  of  c^  with  a  fixed  value  of  C2  =  c  =  0.2.  In  a 
graph  of  N^,  i  =  1  and  2,  versus  p^ ,  Figure  6,  finite  element  results  are 
shown  as  dots.  It  can  be  seen  that  the  dots  may  be  joined  approximately  by 
straight  lines.  Similar  results  are  also  obtained  and  shown  in  the  same  graph 
for  C2  *  c  =  0.3. 

For  three  radial  cracks,  the  crack  2  and  crack  3  remain  the  mirror  image 
of  each  other  so  that  the  problem  is  symmetric  with  respect  to  the  x-axis. 

Two  subcases  are  considered.  In  the  first  case,  the  crack  depths  of  both 
crack  2  and  crack  3  are  fixed  at  C2  =  c]  =  c  =  0.2  or  0.3.  Stress  intensity 
factors  at  crack  tips  of  crack  1  and  crack  2  are  computed  for  various  values 
of  pi*  Finite  element  results  are  shown  as  dots  for  c  =  0.2  and  as  crosses 
for  c  =  0.3  in  Figure  7.  In  the  second  case,  the  first  crack  has  a  fixed 
crack  depth  and  the  other  two  cracks  grow  simultaneously  at  the  same  rate. 
Finite  element  results  for  the  latter  case  arc  used  only  for  the  purpose  of 
checking  our  numerical  approximations. 

In  the  case  of  four  radial  cracks,  crack  U  is  kept  as  the  mirror  image  of 
crack  2.  This  symmetry  reduces  the  problem  to  the  upper  half  of  the  ring  with 
three  cracks.  There  are  three  subcases  which  have  been  studied.  The  first 


case  is  that  only  crack  1  grows.  In  the  second  subcase,  crack  1  and  crack  3 


are  growing  simultaneously  at  the  same  rate  and  crack  2  has  a  fixed  crack 
depth.  The  third  case  is  just  opposite  to  the  first  case.  The  crack  1  is 
fixed  in  length  with  crack  2  and  crack  3  growing  at  the  same  rate.  Using  the 
idealization  shown  in  Figure  4,  finite  element  results  of  stress  intensity 
factors  are  obtained  at  each  crack  tip  for  various  increments  in  crack  depth 
of  the  growing  crack  or  cracks.  The  results  for  the  first  subcase  are  shown 
as  dots  and  crosses  in  Figure  8.  The  numerical  results  for  the  other  two 
subcases  are  not  given  since  they  are  used  only  for  checking  the  numerical 
approximations  by  total  differentials. 


APPROXIMATION  BY  TOTAL  DIFFERENTIALS 

We  assume  that  the  differential  change  in  stress  Intensity  factor  at 
crack  i  due  to  the  differential  change  in  length  of  crack  j  may  be  given  by 

3Nj_ 

dN|  * - dp <  (2) 

3pj 

The  change  in  stress  intensity  factor  due  to  a  small  change  from  pj  to  pj  + 

Apj  can  be  expressed  by 


Pj+APj  3Ni 

ANi  *  /  dpj 

Pj  3pj 


(3) 


If  the  Increment  Apj  is  divided  into  mj  intervals 

Apj  =  A^pj  +  A2Pj  +  ...  +  AmjPj  (4) 

and  if  in  any  interval  Z  there  exists  a  quantity  (3Nf/9pj)£  such  that 


3Nt  Pj+A£Pj  3Nt 

Z  AfcPj  =  /  dpj 

3pj  Pj+Ajt_].pj  3pj 


(5) 
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The  Integral  at  the  right  hand  side  of  Eq .  (3)  may  he  replaced  by  a  summation 

m1  3Ni 

ANl  =  I  (~ — )  (6) 

£=1  3pj 

Furthermore,  If  the  change  In  stress  Intensity  factor  at  crack  tip  1  due  to 
simultaneous  changes  In  length  of  two  or  more  cracks  may  be  approximated  as 
the  sum  of  changes  due  to  each  c rack  length  change,  then  the  total  change  of 
stress  intensity  factor  at  crack  i  due  to  length  changes  of  all  cracks  may  be 
approximately  expressed  by 

n  3Nt 

=  \  1  (“ — )  (7) 

j=l  1=1  3pj 

The  partial  derivative  (3N^/3pj){,  may  be  calculated  from  the  slope  of  the 
tangent  to  the  curve  of  versus  Pj  which  is  obtained  by  curve  fitting  of 

discreted  points  in  the  p-N  plane  from  the  finite  element  computations.  The 

change  in  Mi  must  be  restricted  to  that  due  to  the  change  in  pj  only.  Figures 
6  through  8  are  used  for  this  purpose  with  finite  element  values  of  Mi 

obtained  for  various  values  of  Pi  with  the  crack  depth  of  all  other  cracks 

fixed.  A  study  of  data  points  (p,N)  in  Figures  6  through  8  shows  that  a 
straight  line  of  the  form 

N-l  =  ra(p-l)  (8) 

fits  these  points  in  a  close  approximation.  The  slope  m  of  a  least  square 
regression  Is  given  by 

m  =  I(Ni-l)(Pi-l)/V(Pi-l)2  (9) 

i  i 

This  slope  may  be  used  as  the  average  value  of  the  partial  derivatives  SNj/Spi 

in  a  proper  range  of  pi  for  a  given  set  of  values  n  and  c.  For  n  =  2,  3,  and 

4,  the  slopes  m  of  linear  regressions  are  computed  and  shown  in  Figures  6 


7 


through  8  for  c  =  0.2  or  c  -  0.3  in  a  range  1  «•  ^  1«8 

For  a  given  equal  depth  crack  configuration  we  have 


DNj^ 

3N 

II 

t  — ' 

1  "1 

1 

Dp  ^ 

Dp, 

Dp 

( i ,  j  =  1.2,.. .n) 


(10) 


Stress  intensity  factor  at  any  crack  tip  of  unequal  radial  cracks  may  be 
estimated  by  using  Eqs.  (7)  and  (10).  The  approximate  method  is  applied  to 
subcase  2  of  n  =  3  and  subcases  2  and  3  of  0=4.  The  results  are  compared  to 
finite  element  results  and  the  difference  is  in  general  less  than  four 
percent.  Two  examples  of  the  comparison  are  given  in  Table  I.  In  both 
examples  the  depth  of  crack  1  was  held  constant  and  the  depth  of  all  other 
cracks  increased  an  amount  of  forty  percent  of  the  depth  of  crack  1. 


TABLE  I.  COMPARISON  OF  Kj/p/R-f,  FINITE  ELEMENT  RESULTS  VERSUS  APPROXIMATIONS 


No.  of  Crack 

Crack  Denth 

i  =  i 

:T:r:; 

i 

=  3 

n 

ci 

Other  c 

F.E. 

Anprox. 

BB 

F.E. 

Approx. 

3 

0.2 

0 . 28 

2.26 

2.29 

2.69 

2.70 

o-. 

* 

_ 

2.70 

4 

0.3 

0.42 

2.32 

_ 

2.32 

3.60 

3.50 

3.02 

2.97 

In  the  previous  example  of  four-crack  case,  the  stress  intensity  factor 
of  3.6  at  crack  2  or  4  is  well  above  3.3,  the  would  be  value  If  crack  l  has 
not  arrested  but  grew  at  the  same  rate.  The  stress  Intensity  factor  at  crack 
2  due  to  the  presence  of  crack  1  and  3  Is  considerably  lower  than  4.2,  the 
stress  intensity  factor  for  two-crack  case  with  c  =  0.42.  The  faster  Increase 
In  stress  intensity  factor  at  crack  2  and  crack  4,  however,  tends  to  change 
the  four-crack  configuration  to  a  crack  configuration  dominated  by  two  cracks. 
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The  approximate  method  may  be  used  to  obtain  results  for  crank 

configurations  which  are  more  difficult  ro  compute  by  the  finite  element 

method.  For  instance,  a  ring  with  three  radial  cracks  of  different  depths, 
say  Cj_  ■=  0.2,  C2  =  0.26,  and  cq  =  0.3,  the  entire  ring  with  three  crack  tips 

has  to  he  considered  by  the  finite  element  method.  The  method  of  total 

differentials  can  easily  yield  the  following  estimates:  =  1.026,  - 

1.166,  and  Nq  =  1.26.  Taking  Ke/p/Rq~  =  2.23  from  Figure  6  for  n  =  3,  c  =  0.2, 
we  have  Kq/p/lVq  -  2.28,  ^/p^Ri  -  2.60,  and  Kq/p/Rq  =  2.81. 

For  a  given  value  of  n,  if  we  obtain  PN/Jp  for  many  values  of  c,  then  a 

plot  of  vN/bp  versus  c.  may  be  obtained.  Such  a  plot  for  n  -  3  is  shown  in 

Figure  9.  It  shows  that  JN^/3pj  is  nearly  a  linear  function  of  c.  Let 

be  a  quadratic  function  of  c  of  the  form 

tii^/^Wl  =  :n'c2  (11) 

The  curve  which  best  fits  data  in  the  least  square  sense  is  given  by 

m'  =  c12(c»N2/dpL)i/^.c14  (12) 

i  i 

For  the  curve  in  Figure  9,  m'  is  approximately  0.577.  The  equation  ot  the 
linear  regression  in  Figure  9  is 

dNq/dpq  -  1.077  c  +■  0.29h  (13; 

A  plot  such  as  Figure  9  or  equations  similar  to  Eqs.  (li)  and  (13)  may  be  used 
to  obtain  approximate  values  of  oNj_/ dpi  at  a  specific  value  of  c  for  which  no 
finite  element  results  are  available.  Using  SNj/Cvi  =  0.4037  and  vl^/^Pi  = 
5.77  x  10  for  n  =  3,  c  =  0.1  in  Eqs.  (7)  and  (10)  and  taking  Kg/pr^Y  =  1 • 5h 
from  Figure  5,  our  estimations  are  Kq  /p/Rf  =  1.563,  =  1.438  and 

K3/p^T  =  1  *873  for  the  three-crack  case  of  cj  =  0.1,  c2  =  0.08,  and  cj  = 
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0.15.  In  our  computation  a  negative  value  of  -0.2  was  used  for  Ap2  In  Eq .  (7) 
to  Indicate  the  decrease  In  length  of  crack  2.  The  final  estimation  of  stress 
Intensity  factors  depends  slightly  on  the  judicious  choice  of  c.  In  general, 
we  take  the  smallest  value  of  all  crack  lengths  as  c.  In  the  previous 
example.  If  c  3  0,08  Is  used,  then  c.3  would  have  a  length  Increment  of  more 
than  80  percent.  The  estimations  would  have  been  K^/p/R^  =  1.538,  K^/p^Rq  = 
1.406,  and  K3/ p/ R3  =  1.869.  The  deviations  In  this  example  are  relatively 
small.  In  some  other  cases  the  deviations  become  quite  large  when  a  length 
change  Is  excessive.  The  choice  of  c  for  shallow  cracks  Is  more  Important  and 
will  be  discussed  in  the  following  section. 

APPROXIMATIONS  FOR  SHALLOW  CRACKS 

The  finite  element  method  has  some  difficulty  In  calculating  stress 
Intensity  factors  for  very  shallow  cracks  due  to  the  small  size  of  crack  tip 
elements  required  to  solve  such  problems.  For  c  -  0.05,  finite  element 
results  are  available  In  Reference  16  for  various  n.  For  shallow  cracks  with 
c  <  0.05,  no  finite  element  results  are  available  since  the  accuracy  of  such 
computations  Is  questionable.  Some  results  for  c  <  0.05  were  obtained  by 
other  methods.  The  most  accurate  results  were  reported  In  Reference  20  by 
MMC,  the  modified  mapping  collocation  method.  However,  the  approximate 
methods  for  shallow  cracks  discussed  In  Reference  21  are  easy  to  use  and  yield 
reasonably  accurate  results.  The  crosses  In  the  range  0  <  c  <  0.05  In  Figure 
5  are  obtained  from  the  approximate  formula  (ref  22)  for  a  single  crack 

K  =  {1.12  ao(r=Rj)  “  0 .68(  o0(r=R3 )  -  o0(r=rc)  ]  }/itc  (14) 
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where  a(J(r=Rf)  and  Ojj(r=rc)  are  hoop  stressor  of  an  uuc racked  cylinder  at  r  - 
and  r  =  r  . ,  respectively.  rc  denotes  the  radius  of  the  crack  tip. 

Using  them  and  finite  element  results  for  c  -  0.05  as  a  guide,  the  curves  in 
Figure  5  are  extended  In  broken  lines  to  the  region  c  <  0.1.  The  extended 
curves  agree  well  with  available  'IMG  results. 

Using  the  approximate  method  described  in  the  preceding  section  and 
values  from  Figures  5  and  9,  stress  intensity  factors  can  be  easily  estimated 
for  equally  spaced  shallow  cracks.  The  judicious  choice  of  c  in  shallow  crack 
cases  is  very  important.  For  small  o  and  small  c  the  crack  Interaction  is 
weak.  it  is  recommended  to  choose  c  -  cj  for  the  computation  of  Kj .  For 
instance,  n  =  3,  C]_  •=  0.01,  c2  =  0.02,  and  eg  =  0.03,  we  obtain  Kj/p.lTj  -  0.52 
by  assuming  the  final  crack  configuration  is  reached  from  C|  =  C2  =  C3  -  c  = 
0.91.  In  tiie  computation  of  K2 >  we  assume  the  final  crack  configuration  is 
reached  from  cj  =  c2  =03  =  0-  0.02.  The  result  is  Ki/p/Rj  =  0.73.  For  Kg, 
initial  crack  length  c  -  0.03  is  assumed,  depths  of  crick  1  and  crack  2  are 
decreased  from  0.03  to  eg  ■=  0.01  and  C2  =  0.02,  the  computed  result  is 
K.3/ p/ R'{'  =  0,00.  The  corresponding  results  by  using  Kq .  (Ik)  are  Kq/p/Rg  = 
0.53,  :<2/p/Rf  =  0.7k,  and  Kg/p/R]  =  0.90.  Another  easier  alternative  mpthod 
for  shallow  cracks  is  the  use  of  stress  intensity  f  tutors  for  equal-depth 
cracks.  For  n  =  3,  when  all  crack  depths  are  c  ■=  0.!'l,  the  stress  intensity 
factor  at  a  crack  tin  Is  K/p/Rj  =  0.52.  Similarly  K/p>'Rj  =  0.73  for  c  -  0.02 
and  K/p/Rf  =  0.90  for  c  -  0.03.  Using  these  values  to  the  crick  of  the  right 
crack  depth  in  the  unequal-depth  crack  conf Igurat ion,  all  three  approximate 
methods  for  shallow  cracks  are  useful  up  to  e  =  0.1.  As  a  last  example,  n  = 

3,  eg  =  0.05,  c2  =  0.07,  and  eg  =  0.09,  the  first  method  gives  Kg/p/iTi  ~ 


1  1 


1.122,  K2 / P"^ R-L  "  1«31,  and  K3/p/Rj~=  1.495.  Equation  (14)  gives  K^/p/Rj"  = 
1.15,  K2/P^T  =  1-347  ,  and  K3/p/R'f  =  1.512.  The  use  of  equal  crack  depth 
results  gives  K^/p/R^"  =  1.12,  Ka/p^R^  =  1.31,  and  K3/p/R3  =  1.50. 


CONCLUSIONS 

The  12-node  quadrilateral  Isoparametric  elements  with  collapsed  singular 
elements  around  a  crack  tip  can  be  used  to  efficiently  compute  stress 
intensity  factors  at  tips  of  multiple  radial  cracks  of  unequal  depths  in  a 
thick-walled  cylinder.  Rased  on  finite  element  results  of  some  selected  crack 
configurations,  the  method  of  total  differentials  can  estimate  stress 
intensity  factors  accurately.  The  approximate  method  requires  only  a  small 
number  of  finite  element  computations  to  obtain  graphs  such  as  Figure  7  and 
Figure  9  for  a  given  n,  then  the  method  can  estimate  stress  intensity  factors 
far  an  arbitrary  set  of  unequal  cracks  without  the  need  of  finite  element 
computation  for  that  particular  crack  configuration.  The  advantage  of  the 
approxinate  method  is  not  only  time  saving,  but  also  is  able  to  give  an 
estimate  when  the  problem  is  not  easily  computed  by  the  finite  element 
method . 

For  shallow  cracks,  the  stress  Intensity  factor  at  a  crack  tip  may  be 
estimated  from  equal-depth  crack  configuration  of  that  particular  depth.  This 
is  due  to  the  fact  that  the  change  in  stress  intensity  factor  at  a  shallow 
crack  due  to  changes  in  crack  depths  of  other  shallow  cracks  Is  very  small. 
Alternatively,  each  crick  may  be  considered  as  a  single  crack  since  the 
Interaction  among  shallow  cracks  is  very  weak.  Hence  the  stress  intensity 
factor  at  a  crack  tip  Is  approximately  proportional  to  the  square  root  of  the 
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crack  depth.  A  deeper  crack  will  have  a  higher  stress  intensity  factor  which 
in  turn  will  cause  further  growth  of  the  deeper  crack.  This  process  will 
change  a  large  number  of  shallow  cracks  Into  a  crack  configuration  consisting 
of  several  dominant  cracks.  Further,  if  one  or  two  cracks  are  growing  more 
than  the  rest,  the  increase  in  stress  intensity  factor  at  these  deeper  cracks 
is  greater  than  the  increase  of  others.  This  will  reduce  the  number  of 
dominant  cracks  and  a  single  crack  will  eventually  dominate.  However,  failure 
may  occur  before  reaching  the  final  stage  of  single  crack  domination. 

Although  no  autofrettage  residual  stress  has  been  considered  in  this 
study,  the  approximate  method  aided  by  the  method  summarized  in  Reference  21 
may  be  used  for  partially  autof rettaged,  pressurized,  multiply-cracked 
cylinders  when  crack  depths  are  unequal. 
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10  NODES  COINCIDE 
AT  A  CRACK  TIP 


Figure  2.  Finite  element  idealization  for  two  diametrically  opposite  cracks 
of  unequal  depths. 
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Figure  3.  Finite  element  idealization  for  a  symmetric,  three-crack  problem. 
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DETAIL  B 


DETAIL  A 


Figure  4.  Finite  element  idealization  for  a  symmetric,  four-crack  problem. 
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ATTN:  TECH  LIB  -  DRSTA-TSL 

WARREN,  MI  48090 

COMMANDER 

US  ARMY  TANK-AUTMV  COMD  1 

ATTN:  DRSTA-RC 

WARREN,  MI  48090 

COMMANDER 

US  MILITARY  ACADEMY 

ATTN:  CHMN,  MECH  ENGR  DEPT  1 

WEST  POINT,  NY  10996 

US  ARMY  MISSILE  COMD 
REDSTONE  SCIENTIFIC  INFO  CTR  2 

ATTN:  DOCUMENTS  SECT,  BLDG.  4484 

REDSTONE  ARSENAL,  AL  35898 

COMMANDER 

US  ARMY  FGN  SCIENCE  &  TECH  CTR 
ATTN:  DRXST-SD  1 

220  7TH  STREET,  N.E. 

CHARLOTTESVILLE,  VA  22901 


NOTE:  PLEASE  NOTIFY  COMMANDER,  ARMAMENT  RESEARCH  AND  DEVELOPMENT  CENTER, 
US  ARMY  AMCCOM,  ATTN;  BENET  WEAPONS  LABORATORY,  DRSMC-LCB-TL, 
WATERVLIET,  NY  12189,  OF  ANY  ADDRESS  CHANGES. 
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NO.  OF 
COPIES 


COMMANDER 

US  ARMY  MATERIALS  &  MECHANICS 
RESEARCH  CENTER 
ATTN:  TECH  LIB  -  DRXMR-PL 
WATERTOWN,  MA  01272 

COMMANDER 

US  ARMY  RESEARCH  OFFICE 
ATTN:  CHIEF,  IPO 
P.O.  BOX  12211 

RESEARCH  TRIANGLE  PARK,  NC  27709 
COMMANDER 

US  ARMY  HARRY  DIAMOND  LAB 
ATTN:  TECH  LIB 
2800  POWDER  MILL  ROAD 
ADELPHIA,  MD  20783 


DIRECTOR 

US  NAVAL  RESEARCH  LAB 
2  ATTN:  DIR,  MECH  DIV 

CODE  26-27,  (DOC  LIB) 
WASHINGTON,  D.C.  20375 

COMMANDER 

AIR  FORCE  ARMAMENT  LABORATORY 
1  ATTN:  AFATL/DU 

afatl/dljg 

EGLIN  AFB,  FL  32542 

METALS  &  CERAMICS  INFO  CTR 
BATTELLE  COLUMBUS  LAB 
1  505  KING  AVENUE 

COLUMBUS,  OH  43201 


COMMANDER 

NAVAL  SURFACE  WEAPONS  CTR 
ATTN:  TECHNICAL  LIBRARY  l 

CODE  X212 

DAHLGREN,  VA  22448 


NOTE:  PLEASE  NOTIFY  COMMANDER,  ARMAMENT  RESEARCH  AND  DEVELOPMENT  CENTER, 
US  ARMY  AMCCOM,  ATTN:  BENET  WEAPONS  LABORATORY,  DRSMC-LCB-TL, 
WATERVLIET,  NY  12189,  OF  ANY  ADDRESS  CHANGES. 


NO.  OF 
COPIES 


l 

1 


l 

1 


1 


